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q ■ Abstract 



In a previous study [5] we investigate the bound states of the Hamiltonian de- 
scribing a quantum particle living on three dimensional straight strip of width d. 



We impose the Neumann boundary condition on a disc window of radius a and 



Dirichlet boundary conditions on the remained part of the boundary of the strip. 
We proved that such system exhibits discrete eigenvalues below the essential spec- 
trum for any a > 0. In the present work we study the effect of a magnetic filed of 
Aharonov-Bohm type when the magnetic field is turned on this system. Precisely 
we prove that in the presence of such magnetic filed there is some critical values of 
ao > 0, for which we have absence of the discrete spectrum for < — < qq. We 



d 

give a sufficient condition for the existence of discrete eigenvalues. 
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1 Itroduction 

The task of finding eigenenergies E n and corresponding eigenfunctions / n (r), n = 1,2, ... 
of the Laplacian in the two- (2D) and three-dimensional (3D) domain Q with mixed 
Dirichlet 

/nWlflO D =0 (1-1) 

and Neumann 

nVAMIao =0, (1-2) 
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boundary conditions on its confining surface (for 3D) or line (for 2D) dQ = dQo U dQjy 
(n is a unit normal vector to dQ) [91 [231 ESI HOI El E3] is commonly referred to as 
Zaremba problem [SB], it is a known mathematical problem science. Apart from the purely 
mathematical interest, an analysis of such solutions is of a large practical significance 
as they describe miscellaneous physical systems. For example, the temperature T of 
the solid ball floating in the icewater obeys the Neumann condition on the part of the 
boundary which is in the air while the underwater section of the body imposes on T 
the Dirichlet demand [23]. Mixed boundary conditions were applied for the study of the 
spectral properties of the quantized barrier billiards and of the ray splitting in a variety 
of physical situations. The problem of the Neumann disc in the Dirichlet plane emerges 
naturally in electrostatics [M]. In the limit of the vanishing Dirichlet part of the border 
the reciprocal of the first eigenvalue describes the mean first passage time of Brownian 
motion to dQr>. In cellular biology, the study of the diffusive motion of ions or molecules 
in neurobiological microstructures essentially employs the combination of these two types 
of the boundary coniditons on the different parts of the confinement [33] . 

One class of Zaremba geometries that recently received a lot of attention from 
mathematicians and physicists are 2D and 3D straight and bent quantum wave guides 
[51 H31 UHl ED ESI S31 ST]. In particular, the conditions for the existence of the 
bound states and resonances in such classically unbound system were considered for the 
miscellaneous permutations of the Dirichlet and Neumann domains [2U [2H1 HI] - Bound 
states lying below the essential spectrum of the corresponding straight part were pre- 
dicted to exist for the curved 2D channel if its inner and outer interfaces support the 
Dirichlet and Neumann requirements, respectively, and not for the opposite configuration 
[22~1 [391 S3]- This was an extension of the previous theoretical studies of the existence of 
the bound states for the pure Dirichlet bent wave guide [27] that were confirmed experi- 
mentally. Also, for the 2D straight Dirichlet wave guide the existence of the bound state 
below the essential spectrum was predicted when the Neumann window is placed on its 
confining surface [161 129] . From practical point of view, such configuration can be realized 
in the form of the two window-coupled semiconductor channels of equal widths [291 [30] 
whose experimental creation and study has been made possible due to the advances of 
the modern growth nanotechnologies. The number of the bound states increases with 
the window length a and their energies are monotonically decreasing functions of a [T3] . 
In particular, for small values of a the eigenvalue emerges from the continuous spectrum 
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proportionally to a 4 . The asymptotical estimate for small a were established in [31]. The 
asymptotics expansion of the emerging eigenvalue for small a was constructed formally in 
|5U] , while the rigorous results were obtained in [32] • Recently, this result was extended 
to the case of the 3D spatial Dirichlet duct with circular Neumann disc [5] for which a 
proof of the bound state existence was confirmed, the number of discrete eigenvalues as a 
function of the disc radius a was evaluated and their asymptotics for the large a was given. 
As mentioned above, such Zaremba configuration is indispensable for the investigation of 
the electrostatic phenomena [M]. Similar to the 2D case, it can be also considered as the 
equal widths limit of the two 3D coupled Dirichlet ducts of, in general, different widths 
with the window in their common boundary [7J |3U] . Another motivation stems from the 
phenomenological Ginzburg-Landau theory of superconductivity [2U] which states that 
the boundary condition for the order parameter \I/(r) of the superconducting electrons 
reads. 

The study of quantum waves on quantum waveguide has gained much interest and has 
been intensively studied during the last years for their important physical consequences. 
The main reason is that they represent an interesting physical effect with important 
applications in nanophysical devices, but also in flat electromagnetic waveguide. 
Exner et al. have done seminal works in this field. They obtained results in different 
contexts, we quote [2U [271 ESI EH]- Also in J37J [3SJ H5] research has been conducted in 
this area; the first is about the discrete case and the two others for deals with the random 
quantum waveguide. 

It should be noticed that the spectral properties essentially depends on the geometry 
of the waveguide, in particular, the existence of a bound states induced by curvature 
[T6| I2TI [2U [27] or by coupling of straight waveguides through windows [27] were shown. 

On the other hand, the results on the discrete spectrum of a magnetic Schrodinger 
operator in waveguide-type domains are scarce. A planar quantum waveguide with con- 
stant magnetic field and a potential well is studied in [25] , where it was proved that if the 
potential well is purely attractive, then at least one bound state will appear for any value 
of the magnetic field. Stability of the bottom of the spectrum of a magnetic Schrodinger 
operator was also studied in [57]. Magnetic field influence on the Dirichlet-Neumann 
structures was analyzed in [12 HI], the first dealing with a smooth compactly supported 
field as well as with the Aharonov-Bohm field in a two dimensional strip and second with 
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perpendicular homogeneous magnetic field in the quasi dimensional. 

Despite numerous investigations of quantum waveguides during last few years, many 
questions remain to be answered, this concerns, in particulier, effects of external fields, 
most attention has been paid to magnitec fields, either perpendicular to the waveguides 
plane or threaded through the tube, while the influence of the an Aharonov-Bohm mag- 
netic field alone remained mostly untreated. 

In their celebrated 1959 paper [4] Aharonov and Bohm pointed out that while the 
fundamental equations of motion in classical mechanics can always be expressed in terms 
of field alone, in quantum mechanics the canonical formalism is necessary, and as a result, 
the potentials cannot be eliminated from the basic equations. They proposed several 
experiments and showed that an electron can be influenced by the potentials even if no 
fields acts upon it. More precisely, in a field- free multiply-connected region of space, 
the physical properties of a system depend on the potentials through the gauge-invariant 
quantity § Adl, where A represents the vector potential. Moreover, the Aharonov-Bohm 
effect only exists in the multiply-connected region of space. The Aharonov-Bohm 
experiment allows in principle to measure the decomposition into homotopy classes of the 
quantum mechanical propagator. 

A possible next generalization are waveguides with combined Dirichlet and Neumann 
boundary conditions on different parts of the boundary with a Aharonov-Bohm magnetic 
field with the flux 27ra. The presence of different boundary conditions and Aharonov- 
Bohm magnetic field also gives rise to nontrivial properties like the existence of bound 
states this question is the main objet of the paper. The rest of the paper is organized as 
follows, in Section 2, we define the model and recall some known results. In section 3, 
we present the main result of this note followed by a discussion. Section 4 is devoted for 
numerical computations. 

2 The model 

The system we are going to study is given in Fig 1. We consider a Schrodinger particle 
whose motion is confined to a pair of parallel plans of width d. For simplicity, we assume 
that they are placed at z = and z = d. We shall denote this configuration space by Q 

n = M 2 x [0,d]. 
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Dirichlet boundary condtion 




Figure 1: The waveguide with a disc window and two different boundaries conditions. 

Let 7(a) be a disc of radius a, without loss of generality we assume that the center of 
7(a) is the point (0,0,0); 

7(a) = {(x, y, 0) G M 3 ; x 2 + y 2 < a 2 }. (2.3) 

We set T = dflo\l{ a )- We consider Dirichlet boundary condition on T and Neumann 
boundary condition in 7(a). 
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2.1 The Hamiltonian 

Let be the multiply-connected region fl = {(x,y,z) G fio; x 2 + y 2 > 0}. Let us define, 
now the self-adjoint operator on L 2 (Q) corresponding to the particle Hamiltonian H. This 
is will be done by the mean of Friedrichs extension theorem. Precisely, let Hab be The 
Aharonov-Bohm Schrodinger operator in L 2 (f2), defined initially on the domain C^°(f2), 
and given by the expression 

H A b = (iV + A) 2 , (2.4) 

where A is a magnetic vector potential for the Aharonov-Bohm magnetic field B, and 
given by 

A(x, y, z) = (A 1 , A 2 , A 3 ) = a ( -5-7-5, -f^, o\ a e K \ Z. (2.5) 

\x z + y z x z + y z J 

The magnetic field B : R 3 — > M 3 given by 

B(x, y, z) = curlA = 

outside the z-axis and 

A = 2na, (2.6) 

where 7 is a properly oriented closed curve which encloses the z-axis. It can be shown that 
Hab has a four-parameter family of self-adjoint extensions which is constructed by means 
of von Neumanns extension theory [15]. Here we are only interested in the Friedrichs 
extension of Hab on L 2 (f2) which we now construct by means of quadratic forms. 
For A = (Ai, A 2 , 0) in (j23jl . we observe that A X ,A 2 E L^ c (fi). Let 

tt n = 5(0, n) x [0, d] \ (B(0, l/n) x [0, d]) , n > 2, 

where B(0,r) denotes the disk with centre and radius r. We define on L 2 (f2 n ) (for each 
n > 2) the quadratic form 

q " [uA = L ( f I + Aiu ) ( j + ■ 4i ") + 1 ('I + 01 + j42 ") 

on the domain 

Q(g n ) = {nGH 1 (fi n ); n[r n = 0}, 
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where H 1 (f2 n ) = {u G L 2 (f2 n )|Vw G L 2 (f2 n )} is the standard Sobolev space, T n = dQ n \ 
7(a) and we denote by u\T n , the trace of the function u on T n . The form q n is closed, 
since A l7 A 2 G L°°(ft n ). 
Let the quadratic form g 

q[u,v] = q n [u,v] if «, v G Q(g n ), 
Q(g) = u n ,Q(g n ). 



Lemma 2.1 TTie /orm g closable. 

Proof. The form q is closable if and only if any sequence, u n G Q(q), such that 

lim || u n ||= and lim q[u n — u m ] = 0, (2-7) 



m,n—> 00 



satisfies lim^oo g[w n ] = 0. 
Observe that (12.71) implies 

C := supg[n„] 1/2 < 00. (2.8) 

n 

Let ba a sequence, (n n ) n G Q(g), such that (12. 7p satisfies, then we take e > and choose 
Uq such that 

q[u n — u m ] < e for n,m>n , (2.9) 

and 

II u n \\< e for n > Uq. (2-10) 

Set, moreover, K = fl no C O such that the support of u no included in K. In view of (12.71) 
it follows that 

/ I (z'V + A)(u n — Um) | 2 < q[u n — u m ] — > as n,m— >■ 00, (2-H) 
Jk 

u n | 2 — > as n -> 00. (2-12) 



A is bounded on as 

II A || < I if I n , (2.13) 

we obtain that 

/ I Au n \ 2 — > as n ->• 00. (2.14) 
Since the norm in L 2 is 1-Lipschitz, then 
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1/2 / r \ 1/2 

Mu n - u m ) P dx I - I / I V(u n - u m ) I 2 cfe 



' K J \JK 

<(J | (iV + A)(u n - u m ) | 2 cfe^ . 

According to (I2.12p . the first term on the left side of the latter tends to zero as n, m — > oo 
and, due to ( 12. 9p . the same holds for the right side. Thus, 

/ | u n — u m | 2 + | V(w n — u m ) | 2 — > as n, m — > oo. (2-15) 

Since the form of the classical Dirichlet-Neumann Laplacian in Q(q no ) is closable it follows 
from (I2.15p . in conjunction with (12. lip that 

|2 , n / I „ |2 



Let's consider the following quadratic form 

q[u n ] = q[un,u n ~ % ] + u n ] < q[u n ] 1/2 q[u n - w„ ] 1/2 + | g[w n - u no ] \ . (2.17) 

It follows from (IHD and that 

g[u n ] 1/2 g[u n - u no ] 1/2 < Ce 1/2 when n > n . (2.18) 

Since A is bounded on K we infer from (I2.14p and (12.161) that 

q[u n ,u na ]= / (z'V + A)u n (iW + A)u no -> as n oo. (2.19) 

Using (l2.18p -f l2TT9|) in (I2.17P shows that lim^oo q[u n ] = 0, this ends the proof of the 
lemma |2~T1 ■ 

We denote the closure of q by q and the associated semi-bounded self-adjoint operator 
is the Friedrichs extension of Hab is denoted by H and its domain by D(Q). It is the 
hamiltonian describing our system. We conclude that the domain D(Q) of H is 

D(Q) = {neH'fO); (iV + A) 2 m G L 2 (fi),w[r = 0,z/.(iV + A)w[ 7 (a) = 0}, 

where v the normal vector and 

Hu = (iV + A) 2 n, Vu G D(Q). (2.20) 



2.2 Some known facts 



Let's start this subsection by recalling that in the particular case when a = 0, we get 
H°, the magnetic Dirichlet Laplacian, and when a = +00 we get H°°, the magnetic 
Dirichlet-Neumann Laplacian. 

Proposition 2.2 The spectrum of H° is [(5) 2 ,+oo[, and the spectrum of H°° coincides 



7T \2 



with [{ 2dJ , 
Proof. Since 



00 . 



H = {iV + Af® I® I® (-A M ), on L 2 (M 2 \{0})®L 2 ([0,d]). 



where A := a [ x i y +y i , x 2+ y i ) • Consider the quadratic form 



q[u] = / I (iV + A) it | 2 dxdy 



id x + a- 



y 



u 



x 2 + y 2 

By introducing polar coordinates we get 



If 



dxdy + 



x a 

— = cosy, 

r 



x 

id y - a— 

x 2 + y 2 



y ■ a 
— = sin v. 



II 



dxdy. (2.2i; 



and 



dx 



-y 86 x 



d yd 



= — , d x = cos6- - — , d y = sm9— + 



d x d 



dy r 2 



dr r 2 89 



Hence f!2.2ip becomes 
q[u] = 



d yd sin 9 

1 cos 9- ^7^7: + a I u 

or r z o9 r 

d . x d cos 9 

ism9— + ^ttt: - ct I u 

or r z o9 r 



rdrd9 



rdrd9 



cos 2 9 I d r u I 2 4 



sin 2 ^ 



+ 



(idfjU — au) I 2 I rdrd9 



cos 2 9 \ 
sin 2 ^ I d r u | 2 H — | (idgu — au) | 2 I rdrd9 



d r u 1 2 I (idf)U — au) I 2 | rdrd9 



Expanding u into Fourier series with respect to 9 



u(r,9)= jr 



Jk8 



Uk(r, 



k=— 00 



2tt 



(2.22) 
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enables us to rewrite (I2.22p as 

d r u | 2 +— | (idgu — au) \ 2 ) rdrdO 



> -4 



^2 -( k + 



rdrdO 



I r 2 




\k + 


min 

k 


k + 


a | 2 


min 

k 


k + 


a | 2 


min 

k 


k + 


a | 2 


min 


k + 


a | 2 



7, 



k 



ikd 



—z u\r, 



x 2 + y 2 



2vr 
2 rdrdO 



rdrdO 



u(x, y) | 2 dxdy. 



Hence 



(zV + A)u | 2 dxd?/ > min | k + a 

k 



-— — - | | 2 did?/. (2.23) 



Here the form in the right hand side is considered on the function class H 1 (R 2 ), 
obtained by the completion of the class Cq°(IR 2 \{0}). Inequality (12 . 23[) is the Hardy 
inequality in two dimensions with Aharonov-Bohm vector potential [3]. This yields that 
a ((iV + A) 2 ) C [0,+oo[ 



Since cr(— A) = a ess (— A) = [0, +oo[, then there exists a Weyl sequences {h n }^ =1 for 
the operator —A in L 2 (IR 2 ) at A > 0. Construct the functions 



h n if x > n and y > n, 
if not. 



+ A) 2 - AJ ip n || < || (A - \)<p n || + || AV„ II + II AV<^ r 
< ||(A-AK||+^. 



Where c is positive real. 
Therefore, the functions ip n 



is Weyl sequence for (z'V + A) 2 at A > 0, thus 
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[0, +oo[c a ess ((zV + A) 2 ) C o ((z'V + A) 2 ) . 



Then we get that the spectrum of (z'V + A) 2 is [0, +oo[, we know that the spectrum 



of -A° M] and -Ajfo is {( 



d ' ' 



j e N*} and {( 



7T 

Therefore we have the spectrum of H° is [(-r) 2 , 

7T 



2d 



j G N} respectively. 



-oo[. And the spectrum of H°° coincides 



with [(^r,+°o[. ■ 

Consequently, we have 



,7T, 

~d' 



-oo 



C a(H) C 



7T 



-oo 



Using the property that the essential spectra is preserved under compact perturbation, 
we deduce that the essential spectrum of H is 



,7T, 

I'd' 



+00 



2.3 Preliminary: Cylindrical coordinates 

Let's notice that the system has a cylindrical symmetry, therefore, it is natural to consider 
the cylindrical coordinates system (r, 9,z). Indeed, we have that 

L 2 (Q,dxdydz) = L 2 (]0, +oo[x [0, 2vr[x [0, d], rdrdOdz). 

Consequently, a corresponding orthonormal basis in M 3 is given by the three vectors 




cos 9 \ / — sin 8 

sin 8 , ee = I cos^ 
J \ 

We note by Ag, the Aharonov-Bohm magnetic potential vector (12. 5p in cylindrical coor- 
dinates given by 

Ag(r,8,z) = — (— sin0, cos^, 0) = —eg. 

We denote the gradient in cylindrical coordinates by V r ,e,z- While the operator z'V + A 
in cylindrical coordinates is given by 



iV r .g :Z + Ag = i—e r + - I 1— + a 
or r \ 08 



d 

oz 



thus the Aharonov-Bohm Laplacian operator in cylindrical coordinates is given by 



iV r ,g, z + A 6 



1 d d I d 2 d 2 
r dr dr r 2 88 d 2 z 



11 



3 The result 



The main result of this paper is the following: 



Theorem 3.1 Let H be the operator defined on (I2.20p . There exist ao > such that for 

o- d (H) = 0. 



any < — < a , we have 
a 



There exist a x > 0, such that — > a\, we have 

a 

v d (H) ^ 0. 

Remark 3.1 The presence of magnetic field in three dimensional straight strip of width 

d with the Neumann boundary condition on a disc window of radius < — < Oq and 

a 

Dirichlet boundary conditions on the remained part of the boundary, destroys the creation 

of discrete eigenvalues below the essential spectrum, if — > a\, the effect of the magnetic 

a 

field is reduced. 

Proof. As in [5], let's split L 2 (Q,rdrd6dz) as follows, L 2 (Q,rdrdddz) = 
L 2 (n-,rdrd6dz) © L 2 (fi+, rdrdOdz), with 

n~ = {(r,8,z) g [0,a] x [0,2tt[xM}, 

q+ = n\n- 



Therefore 



Ha' N © H+> N <H< H-' D © Ht' D . 



Here we index by D and N depending on the boundary conditions considered on the 
surface r = a. The min-max principle leads to 

°ess{H) = o- ess (H+' N ) = a ess (H^°) = [C^f \ +oo[. (3.24) 

Let us denote by \k(H~' N ), X k (H~' D ) and Xk(H), the k-th eigenvalue of H~' N , H~' D and 
H respectively then, the minimax principle yields the following 

X k (H~' N ) < X k (H) < X k (H-' D ) (3.25) 

and for k > 2 

X k ^(H-' N ) < X k (H) < X k (H-' N ); (3.26) 
12 



h-i(K> D ) < X k (H) < X k (H^ D ). (3.27) 

Thus, if H~' N does not have a discrete spectrum below (f ) 2 , then H do not has as well. 
We mention that its a sufficient condition. 
The eigenvalue equation is given by 

H;> N f(r,6,z)=Ef(r,6,z). (3.28) 

This equation is solved by separating variables and considering f(r, 9, z) = R(r)P(6)Z(z). 

Plugging the last expression in equation f !3.28|) and first separate Z by putting all the z 

Z" 

dependence in one term so that — can only be constant. The constant is taken as 

— k 2 z = ) 2 f° r convenience. Second, we separate the term t^tj-t + a) 2 P 

which has all the 9 dependance. Using the fact that the problem has an axial symmetry 

1 d 

and the solution has to be 2n periodic and single value in 9, we obtain — {i-^ + &) 2 P 

P o9 

should be a constant — (m — a) 2 = —v 2 for m G Z. Finally, we get the following equation 
for R 

R"(r) + ^R'(r) + [E - k 2 z - ^}R(r) = 0. (3.29) 

We notice that the equation ( 13.291) . is the Bessel equation and its solutions could be ex- 
pressed in terms of Bessel functions. More explicit solutions could be given by considering 
boundary conditions. The solution of the equation (13. 29ft is given by R(r) = cJ v (/3r), with 
c G R, (3 2 = E — k 2 and J u is the Bessel function of first kind of order v. We assume 
that R'(a) = 0. So we get the condition a(3 = x' un , with x' un are the positive zeros of the 
Bessel function J' n . Thus H~' N has a sequence of eigenvalues given by 

A ■ - 4- k 2 

Lb 

x' 2 n /(2j + l)vr x 



a 2 V 2d 



As we are interested for discrete eigenvalues which belongs to [(t^) > (^j ) on ^ -V^r 



intervenes. 
If 



7T 



2 



< Ao,^n, (3.30) 

then there H does not have a discrete spectrum. We recall that v 2 = (m — a) 2 and it is 
related to magnetic flux, also recall that x' v n are the positive zeros of the Bessel function 



13 



J' n and W > 0, Vn 6 N*; < x' v n < x' v n+1 ( see [SI]). So, for any eigenvalue of H a ,7V , 



a 2 V2(i7 a 2 V2d/ 
An immediate consequence of the last inequality is that to satisfy ( 13.30]) it is sufficient to 
have 



therefore 



then 



6 \2d) < a? ' 



V3n <,i 
2d < a 



a 2x' vl 
d v^tt' 



Or ( see [T71 SU SHI [54]), we have 



v + a n z/ 1/3 < x' u n , 
where a n = 2 _1 / 3 /3 n and f3 is the n-th positive root of the equation 



For n = 1, we have 



J | fiz 3 / 2 )-J = a[;x 3/2 l=0. 



c := 0.6538 + a < 0.6538 + i/ < x' u v (3.31) 



Then we get that for d and a positives such that — < : 



u ' v/3V 
M#) = 0- 

This ends the proof of the first result of the theorem 13.11 

By the min-max principle and ( I3.27p . we know that if H~ ,D exhibits a discrete spec- 
trum below (^) 2 , then H do as well. 
H~ ,D has a sequence of eigenvalues [SI ISS], given by 

■x»,n\ 2 , f(2j + l)n 



Aj>,n - ( ) + 



2d 



Where x v , n is is the n— th positive zero of Bessel function of order v ( see [S]). As we are 



7T 7T 

interested for discrete eigenvalues which belongs to [(^) 2 ; (^) 2 ) on ^ ^ or -V^n- 
If we is satisfied the condition 

7T\ 2 



\u,n < [2) , (3-32) 
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then H have a discrete spectrum. 

We recall that < x v>n < x U)n+ \ for any v > and any n G N* ( see [SI]). So, for any 
eigenvalue of H~ ,D , 



a 2 V2(i7 a 2 V2d7 
An immediate consequence of the last inequality is that to satisfy ( I3.32p it is sufficient to 
set then 



2x u ^\ a 



Or ( see [TTJ we have 



1\ 2 

n — - ) n 2 + v 2 < x Utn , 



3ir\* n I /37T X z 



For n = 1, we have 



Then we get that for d and a positives such that — > a x : 



d v/3V 
cr d {H) + 0. 



4 Numerical computations 

This section is devoted to some numerical computations. We represent the radius values 
a : a (->• ^j=^ and ai : a ^j^r, where and a; Q) i is the first positive zero of Bessel 
function J' a and J a respectively, a function of the flux magnetic a, which makes it possible 
to exist or not discrete eigenvalues of H. 
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The values of ^ which appears the eigenvalues of H. 




Figure 2: The waveguide with a disc window and two different boundaries conditions. 



In the neighborhood to zeros, we do to zoom the figure [21 when d — > 0, we observe 
that the discrete eigenvalues of H exists, contrary when a — > 0. 
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Figure 3: Zoom of figure |2] in the neighborhood to zeros. 

We analyse a dependence of the first eigenfunction on the flux magnetic a in the whole 
range of their variation. We chose a more natural method: the finite element method. 
Consider a triangulation of Q. We will construct the matrices of mass M and of 
assemblage A associated to the mesh and to the operator. We consider two elements 
u, V of D(Q). Note U, V vector whose i-th coordinate gives the value of u, v at i- 
th point of the mesh. Denote by the basis functions for the finite element element 
Pfc considered, the matrices M and A are given by their coefficients Mjj and Ajj are 
numerical approximations of integrals: 



M 4 



1,3 



A ■ ■ ~ 



[ (*V + A)0^V + A) 
Jn 



3 ' 
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Matrices M and A satisfied then: 

/ uv « t VMU, 
a(u,v) « *FA?7. 

From calculation of the two matrices M and A, we wish to determine the eigenvector of 
the operator H. We can draw a parallel between the initial theoretical problem and the 
numerical problem that we solve: 

• find u £ D{VL), jiiGl such that Wv £ -D(fi), a(u, t> ) = fi < u, v >. 

• find U £ Ft, n £ R such that £ P fc , *FAt/ = yfVMA. 

The approximate problem is to determine the eigenvectors of the matrix A for the matrix 
M, that is to say, to find U £ ¥ k , such as: AU = XMU. 

We represent the first eigenfunction of H for the Neumann radius a=3, d = n and 
a = 0.5. 








Figure 4: the first eigenfunction of H for the Neumann radius a=3, d = it and a = 0.5. 
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We represent the first eigenfunction of H for the Neumann radius a=3, d = it and 
a — 1. 




Figure 5: the first eigenfunction of H for d = it and a = 0.75. 
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